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Abstract 

The two-body decay rate of a weakly decaying particle (such as the kaon) is shown to 
be proportional to the square of a well-defined transition matrix element in finite volume. 
Contrary to the physical amplitude, the latter can be extracted from finite- volume correla- 
tion functions in euclidean space without analytic continuation. The K — > tttt transitions 
and other non-leptonic decays thus become accessible to established numerical techniques 
in lattice QCD. 



1. Introduction 

The computation of the non-leptonic kaon decay rates from first principles, using 
lattice QCD and numerical simulations, meets a number of technical difficulties (see 
ref. [1], for example). Apart from the operator renormalization, which must be con- 
trolled at the non-perturbative level, the central problem is that the computational 
framework is limited to correlation functions in euclidean space and that there is 
apparently no simple relation between the behaviour of these functions at large time 
separations and the desired transition matrix elements [2,3]. 
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This statement (which is often referred to as the Maiani-Testa no-go theorem) 
applies to very large or infinite lattices, where the spectrum of final states is contin- 
uous. One might think that having a finite volume (as is unavoidable when numeri- 
cal simulations are employed) makes it even more difficult to extract the transition 
amplitudes. In the present paper we wish to show that this is actually not so. The 
key observation is that the two-pion energy spectrum is far from being continuous 
when the lattice is only a few fermis wide. Under these conditions, a kaon at rest 
cannot decay into two pions unless one of these energy levels happens to be close to 
its mass. This is the case for certain lattice sizes, and a simple formula then relates 
the square of the corresponding transition amplitude in finite volume to the physical 
decay rate in infinite volume. 

The problem is thus reduced to calculating the required finite-volume transition 
amplitudes. Since the initial and final states are isolated energy eigenstates, these 
matrix elements can in principle be computed using established techniques, such 
as those commonly employed to determine form factors. An additional difficulty 
is that the relevant two-pion states are not the lowest ones in the specified sector. 
Two-particle states in finite volume have, however, previously been studied [4-14] 
and practical methods have been devised to calculate the higher levels. 

To keep the presentation as transparent as possible, we shall consider a simplified 
generic theory with two kinds of spinless particles, referred to as the kaon and the 
pion. Details are given in the next section, and we then first discuss the form of 
the two-pion energy spectrum in finite volume. This is essentially a summary of 
the relevant results of refs. [15-17]. In sect. 4 we define the transition amplitudes 
in finite volume and state the formula that relates them to the corresponding decay 
rates in infinite volume. The following sections contain the proof of this relation and 
a discussion of its application to the physical kaon decays. 



2. Preliminaries 



As announced above, we consider a generic situation where there are two particles, 
the "kaon" and the "pion" , with spin zero and masses such that 

2m w < rriK < 4/71^. (2-1) 

We assume that the symmetries of the theory are such that the kaon is stable in 
the absence of the weak interactions and that the pions scatter purely elastically 
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below the four-pion threshold. The weak interactions, described by a local effective 
lagrangian £ w (x), then allow the kaon to decay into two pions. The corresponding 
transition amplitude is 

T(K^ttit) = (ir Pl ,irp 2 out\£ w (0)\Kp), (2.2) 

with pi, P2 and p the four- momenta of the pions and the kaon. We shall only be 
interested in the physical case where the total momentum p = p\ + p 2 is conserved. 
Lorentz invariance and the kinematical constraints then imply that the transition 
amplitude is independent of the momentum configuration. 

The meson states in eq. (2.2) are normalized according to the standard relativistic 
conventions (appendix A) and their phases are constrained by the LSZ formalism. 
In the case of the pions, for example, one assumes that there exists an interpolating 
hermitian field (p(x) such that 



(0\<p(x)\7rp) = VZ„e-^ x (2.3) 

for some positive constant Z w . If the phase of the kaon states is chosen in the same 
way, the CPT symmetry implies 

T(K -► vrvr) = Ae iS ° (2.4) 

with A real and So the 5-wave scattering phase shift of the outgoing pion state. The 
decay rate is then given by the usual expression 



proportional to the pion momentum k n in the centre-of-mass frame. 



3. Two-pion states in finite volume 



In a spatial box of size Lx Lx L with periodic boundary conditions, the eigenvalues 
of the total momentum operator are integer multiples of 2-ir/L. The energy spectrum 
is also discrete in this situation, with level spacings that can be appreciable. In the 
following we consider the subspace of states with zero total momentum and trivial 
transformation behaviour under cubic rotations and reflections. 
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The energy spectrum of the two-pion states in this sector below the inelastic 
threshold W = 4m^ has been studied in detail in refs. [15-18]. In particular, for the 
lowest energy value the expansion 



W = 2m n - 



m n L 3 



4-7rao 



{ 



l + Cl ^+c 2 ||} + 0(L- 6 ), 



(3.1) 



d = -2.837297 



c 2 = 6.375183, 



(3.2) 



has been obtained, where 



ao = lim 




(3.3) 



is the 5-wave scattering length (here and below the scattering phase is considered 
to be a function of the pion momentum k in the centre-of-mass frame). The higher 
energy values in the elastic region are determined through 



where n = 1, 2, . . . labels the energy levels in increasing order and the angle 4>(q) is 
a known kinematical function (appendix B). Apart from the lowest level, the energy 
spectrum at any given value of L is thus obtained by inserting the solutions k of 
cq. (3.5) in eq. (3.4) f. 

All these results are valid up to terms that vanish exponentially at large L. Box 
sizes a few times larger than the diameter of the pion should be safe from these cor- 
rections. Equation (3.5) moreover assumes that the scattering phases S[ for angular 
momenta I > 4 are small in the elastic region, which is usually the case since 5i is 
proportional to k 2l+1 at low energies. 

For illustration, let us consider QCD with three flavours of quarks, unbroken 
isospin symmetry and quark masses such that the masses of the charged pions and 
kaons coincide with their physical values. In the subspace with isospin 0, the two- 
pion energy spectrum is then given by eqs. (3.1)-(3.5), with 5$ the appropriate pion 

\ Similar formulae have been derived for the spectrum in the subspaces of states with non-zero 
total momentum [19]. The extension of our results to these sectors could give further insight into 
the connection between finite and infinite volume matrix elements and may prove useful in practice. 



W = 2^Jml + k 2 , 



(3.4) 



mr - S (k) = 4>(q) 



kL 
2^ 



(3.5) 
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Fig. 1. Two-pion energy spectrum in QCD below the inelastic threshold, in the 
sector with isospin 0, calculated from eqs. (3.1)-(3.5) with the scattering phase shift 
given by next-to-leading order chiral perturbation theory. The levels shown in this 
plot are all non-degenerate. 

scattering phase. If we insert the phase shift that is obtained at one-loop order of 
chiral perturbation theory [20-22], this yields the curves shown in fig. 1. For any 
other reasonable choice of the scattering phase the plot would look essentially the 
same, because the interaction effects are proportional to 1/L 3 and thus tend to be 
small. Note that the spacing between successive levels is quite large. One is clearly 
very far away from having a continuous spectrum when L < 10 fm. 



4. Kaon decays in finite and infinite volume 

Let us imagine that a state \K) describing a kaon in finite volume with zero mo- 
mentum has been prepared at time xq = 0. In the absence of the weak interactions, 
this is an energy eigenstate (and thus a stationary state) with energy mx- However, 
through the interaction hamiltonian 



the time evolution of the state becomes non-trivial and it starts to mix with the other 




(4.1) 
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eigenstates of the unperturbed hamiltonian. It is straightforward to work this out 
using ordinary time-dependent perturbation theory. For the transition probability 
at time xq = t to any finite- volume two-pion state \tttt) with energy W, the result 



P(K ^ Ttn) = 4\(int\H w \K)\ 2 ^M", io = W-m K , (4.2) 

is then obtained (in this equation the states are assumed to be normalized to unity 
and higher-order weak- interaction effects have been neglected). 

From eq. (4.2) one infers that the transition probabilities tend to be very small 
unless the energy of one of the two-pion final states happens to be close to the kaon 
mass. Recalling fig. 1, it is clear that this will be the case only for certain box sizes 
L. In the following we focus on these special values of L and introduce the associated 
transition matrix element 



M = (Trir\H v \K), (4.3) 

where both states are normalized to unity as before, while their phase will not matter 
and can be chosen arbitrarily. Since W = mx in this case, eq. (4.2) becomes 

P(K -> tttt) = \M\ 2 t 2 (4.4) 

and the kaon will thus have an appreciable probability to decay into the two-pion 
state if one waits long enough (the formula breaks down at very large times, because 
the higher-order terms are then no longer negligible). 

The central result obtained in the present paper is that the finite-volume matrix 
element M is related to the decay rate of the kaon in infinite volume through 

[cf. eqs. (2. 5), (3. 5)]. The relation holds under the same premises as eq. (3.5) and 
the comments made in sect. 3 thus apply here too. Another restriction is that the 
two-pion final state has to be non-degenerate in the specified sector of the unper- 
turbed theory. This condition is satisfied for n < 8 [17], but degeneracies can occur 
at higher level numbers and the formula then ceases to be valid. 

In principle eq. (4.5) allows one to compute the kaon decay rate in infinite volume 
by studying the theory in finite volume. Note that in the course of such a calculation 
it should also be possible to determine the two-pion energy spectrum and thus the 
scattering phase Sq in the elastic region. 
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The proportionality factor in eq. (4.5) essentially accounts for the different nor- 
malizations of the particle states in finite and infinite volume. One can easily check 
this in the free theory, where the pion self- interactions are neglected. In this case 
and for n < 6, the ro-th two-pion energy level passes through vtlk at 

L = ^v^- (4.6) 

Equation (4.5) then assumes the form 

\A\ 2 = -(m K L) 3 \M\ 2 , (4.7) 

v n = number of integer vectors z with z 2 = n, (4-8) 

which is precisely what is derived from the relative normalizations of the plane waves 
in finite and infinite volume that describe the (non-interacting) kaon and pion states 
(sect. 6). 



5. Proof of eq. (4.5) 

The interpretation of the proportionality factor in eq. (4.5) given above also applies 
in the interacting case. This follows from the fact that the transition matrix elements 
probe the S'-wave component of the two-pion wave function near the origin and that 
this component is the same in finite and infinite volume apart from its phase and 
normalization. The latter can be worked out explicitly in the framework of refs. [15- 
17], but the calculation is rather involved and will not be presented here. 

Instead we shall go through a different argument, where one studies the influence 
of the weak interaction on the energy spectrum in finite volume. This can be done 
directly, using ordinary perturbation theory, or one may start from eq. (3.5) and take 
the weak- interaction effects on the scattering phase into account. The combination 
of the results of these calculations then yields eq. (4.5). 

As already mentioned in sect. 2, the kaon is assumed to carry a quantum number 
(alias strangeness) that forbids its decay into pions in the unperturbed theory. Since 
only the strangeness-changing part of the weak interaction lagrangian contributes 
to the kaon transition amplitudes, all other terms may be dropped without loss. 
The matrix elements of the weak hamiltonian H w between states with the same 
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Fig. 2. Kaon resonance contribution to the elastic pion scattering amplitude in the 
s-channel. The diagram appears at second order of the expansion in powers of the 
weak interaction, with the bubbles representing the first-order Kirir vertex function. 



strangeness are then all equal to zero. As a consequence most energy values in finite 
volume are affected by the weak interaction only to second order. 

First order energy shifts do occur, however, if there are degenerate states at lowest 
order that mix under the action of H w . This is the case at the values of L where one 
of the two-pion energy values coincides with the kaon mass, i.e. at the special points 
considered in the preceding section. Degenerate perturbation theory then yields 

W = m K ±\M\ + ... (5.1) 

for the first order change of these energy values (here and below the ellipses denote 
higher-order terms that do not contribute to the final results). 

The energy shifts (5.1) can also be calculated by including the weak corrections 
to the scattering phase on the left-hand side of eq. (3.5). From the above one infers 
that the solutions of eq. (3.5) we are interested in are given by 

k = k n ±Ak + ..., Ak = ^\M\. (5.2) 

Compared to the kaon resonance width (which is of second order in the weak inter- 
action), these values of k are far away from the kaon pole. The weak corrections to 
the pion scattering amplitude in the relevant range of energies are hence small and 
can be safely computed by working out the perturbation expansion in powers of the 
interaction lagrangian. 

One might think that these corrections are all of second or higher order, because 
the interaction is strangeness-changing. The reason this is not so is that the kaon 
propagator in a diagram like the one shown in fig. 2 evaluates to 

lZK =± iZ * , +■•■ (5-3) 



p 2 — m 2 K 2mx \M\ 

at the energies (5.1) and thus reduces the effective order of the term by 1. This 
diagram is in fact the only one that yields a first-order contribution to the scattering 
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amplitude. It can be calculated by noting that the momenta flowing into the three- 
point vertices are all on shell up to higher-order corrections. The vertices are hence 
proportional to the kaon decay amplitude A. Together with eq. (5.3) this leads to 
the result 

a "» (t) = 4 » (t)T 32^pgi + - < mod " < 54 > 

for the scattering phase in the full theory at the point (5.2) (as in the previous 
section, So stands for the phase shift in the unperturbed theory). 

We now replace So in eq. (3.5) by So and expand all terms in powers of the weak 
interaction. The lowest-order terms cancel while, at first order, the equation implies 

This is easily seen to be equivalent to eq. (4.5) after substituting the expression (5.2) 
for Ak and we have thus proved this relation. 



6. Verification of eq. (4.5) in perturbation theory 



In a low-energy effective theory, such as the chiral non-linear cr-model, it is possible 
to obtain an independent check on eq. (4.5) by working out the transition amplitudes 
in finite and infinite volume in perturbation theory. Since this calculation does not 
rely on any of the results presented above, it can provide additional confidence in 
the correctness of the equation. Perturbation theory may also prove helpful when 
considering more complicated situations, where one has several decay channels or 
particles with non-zero spin. In this section we describe how such a calculation goes, 
without giving too many details. 

6.1 Specification of the model 

The two-pion energy spectrum and the proportionality factor in eq. (4.5) depend on 
the final-state interactions only through the phase shift So- All other properties of 
the pion interactions do not matter and to check the equation we may thus consider 
an arbitrary effective meson theory with the correct particle spectrum. 
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For the pion interaction lagrangian the simplest choice is 



AntO) = ir V(x) 4 , 



(6.1) 



where ip(x) denotes the pion field and A the bare coupling. To make the perturbation 
expansion completely well-defined, we introduce a Pauli-Villars cutoff A. At tree 
level the euclidean pion propagator is then given by 



with m the bare mass of the pion (its physical mass is denoted by m n as before). 
The cutoff should be large enough so that ghost particles cannot be produced at 
energies below the four-pion threshold, but in view of the universality of eq. (4.5) 
there is no need to take A to infinity at the end of the calculation. 

As far as the kaon is concerned, the least complicated possibility is to describe it 
by a hermitian free field 9{x) with mass mx and to take 



as the weak-interaction lagrangian. One then first has to expand the transition am- 
plitude (2.2) in powers of A, but we shall not discuss this here since the calculation 
is completely standard. The way to obtain the perturbation expansion of the finite- 
volume matrix element (4.3) may be less obvious, however, and we thus proceed to 
explain this in some detail. 

6.2 Two-pion states 

In finite volume the low-lying two-pion energy eigenstates with zero total momentum 
and trivial transformation behaviour under the cubic group may be labelled by an 
integer n = 0, 1,2. . . such that the associated energies W n increase monotonically 
with n. We denote these states by \inrn) and assume that they have unit norm. 

To lowest order in A, the energy values are determined through the free energy- 
momentum relation and the relative momentum of the pions. Since we only consider 
cubically invariant states, any two momenta that are related to each other by a cubic 
transformation describe the same state. For n < 6 the momenta in the set 




(6.2) 



£ w (x) = \g9(x)ip(xf 



(6.3) 



n„ = {k = 2vrz/L | z G Z 3 , z 2 = n 



} 



(6.4) 
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Fig. 3. Feynman diagrams contributing to the correlation function (6.7). The lines 
represent the free pion propagator in the time-momentum representation (6.8) and 
the filled circles the self-interaction vertex. All external lines end at times xu or yo. 

are all equivalent in this sense. The corresponding state is thus non-degenerate and 
one concludes from this that 

W„ = 2^m 2 + n(2vr/L) 2 + 0(A), < n < 6. (6.5) 

In the following, our attention will be restricted to these levels. 

The corresponding energy eigenstates \mrn) can be created from the vacuum by 
applying the operators 

O n (x )=Yl f L a 3 xd 3 ye^-y^(x ,xMx ,y). (6.6) 

Note that O n (x (> ) couples to all two-pion states in the given sector, since there are 
no quantum numbers that would forbid this. In euclidean space and at large time 
separations xo — yo, its connected two-point function is thus given by 

6 

(O n (x )O n (y )) con = |<O|0„(O)|7T7rZ>| 2 e - w '(*°-«°> + . . . , (6.7) 
;=o 

where the ellipses stand for more rapidly decaying terms. 

The perturbation expansion of the two-pion energy W n and the associated matrix 
element |(0|C n (0)|7T7rn)| may now be obtained by expanding the left-hand side of 
eq. (6.7) in Feynman diagrams in the standard way. If one uses the time- momentum 
representation 

/ d 3 xe-* px MxM0)> = —z (m~A), Wp = Vm 2 + p 2 , (6.8) 

Jo 

for the tree-level pion propagator, the diagrams evaluate to a sum of exponentials. 
The desired expansions can then be read off from the coefficients of the exponential 
factor that corresponds to the n-th level. 
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Y ? 

a b c 

Fig. 4. Diagrams contributing to the correlation function (6.10). The double line 
represents the kaon propagator and the circled cross the weak interaction vertex at 
the origin. All other graphical elements are as in fig. 3. 

To leading order the diagrams a and b in fig. 3 yield the expected expression (6.5) 
for the two-pion energy and 

|(0|e> n (0)|7T7rn)| = Vl^L 3 /W n + 0(\) (6.9) 

for the matrix element [cf. eq. (4.8)]. At the next order in the coupling, there are two 
types of diagrams. Diagram c and three further diagrams of this kind amount to an 
additive renormalization of the pion mass by a term that is independent of L up to 
exponentially small corrections [15]. Such contributions are neglected here and the 
renormalization is thus equivalent to replacing m by m n in the tree- level expressions. 
One is then left with the diagram d, which can be worked out analytically in a few 
lines. 

6.3 Transition matrix element 

The finite- volume transition matrix element (4.3) can be computed by studying the 
euclidean correlation function 

/ d 6 y (O n (x )Av(0)%)> con — 

JO 

6 

£ e - w ' X0+mKV " (0|O„(0)K l\H w \K)(K\6(0)\0) + . . . (6.10) 

1=0 

at large xq and large negative uq. As in the case of the two-pion states, the terms 
we are interested in are found by looking for the appropriate exponential factor. 
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To lowest order diagram a in fig. 4 yields 



\{irirn\H w \K)\ = ^J^—j {1 + 0(A)} . (6.11) 

The pion mass in this expression is renormalized by the tadpole insertions at the 
next order (diagram b and its mirror image). Diagram c, the only other diagram 
at this order, may be evaluated by inserting the time-momentum representation for 
the external and also the internal lines. Apart from various simple terms, one then 
ends up with the momentum sum 

S n = L~ 3 V { -. I -R A (p 2 ,k 2 )\, ke!l n , (6.12) 

where R\ is an expression that arises from the Pauli-Villars regularization. 

A general summation formula proved in ref. [16] allows one to compute such sums 
up to terms that vanish more rapidly than any power of 1/L. The precise form of 
Ra is not important for this. One only needs to know that it is a smooth function 
of p and k and that it makes the sum absolutely convergent. The result 



d 3 p ( 1 



(2vr) 3 \2u>- 



1 1 

+ 



p 2 — k 2 + it p 2 — k 2 — ie 



flA(p 2 ,k 2 ) 



is then obtained, with the constant z n given by 

z n = lim S^Z 00 (l;q 2 ) + ^-\ (6.14) 
q 2 ^n [ q — n ) 

(the zeta function Z 00 (s; q 2 ) is defined in appendix B). 

6.4 Final steps 

To check eq. (4.5) one has to tune the box size so that W n = tuk for a specified 
level number n. This condition determines L order by order in the coupling. The 
perturbation expansion of the right-hand side of eq. (4.5) is then obtained by insert- 
ing this series in the proportionality factor and the perturbative expressions for the 
matrix element |(7r7rrc|.ff w |.K")|. 
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To lowest order, the box size is given by eq. (4.6) and the function (j)'(q) in the 
proportionality factor is thus to be expanded around q = y/n. This generates a term 
proportional to z n , which cancels the corresponding term in eq. (6.13). The integral 
in this equation matches with the contribution to the transition amplitude A of the 
infinite-volume diagram with the topology of diagram c. All other terms that occur 
at first order in the coupling cancel and one finds that eq. (4.5) holds as expected. 



7. Application to the physical kaon decays 

Compared to the generic theory considered so far, the situation in the case of the 
physical kaon decays is complicated by the fact that there are several decay channels. 
To a first approximation we may however assume that isospin is an exact symmetry 
in the absence of the weak interactions. The decay channels can then be separated 
from each other by passing to a basis of states with definite quantum numbers. 

As an example we discuss the CP-conserving decays of the neutral kaon into two- 
pion states with isospin and 2. The corresponding decay amplitudes, A and A 2 , 
are related to the physical transition matrix elements through 

T(K° S - vr+vr-) = -j=Ao<& + 

T(K° S - vr 7r°) = - -1 A e< + ^= A 2 e iS o . (7.2) 
y6 y3 

In these equations 5$ denotes the 5-wave pion scattering phase in the channel with 
isospin I and the normalization and phase conventions are as in sect. 2. 

In the sector of two-pion states with isospin /, zero electric charge, zero total mo- 
mentum and trivial transformation behaviour under cubic rotations and reflections, 
the energy spectrum in finite volume is determined by the equations that we have 
previously discussed, with 5q replaced by Sq. At the points where one of these energy 
levels passes through m K , we define the associated transition matrix element 

Mj = ((^r) 7 |F w |K°>, (7.3) 

where it is understood that the states are normalized to unity and that H w is the 
CP-conserving part of the effective weak hamiltonian. With these conventions, the 



14 



Table 1. Calculation of the proportionality factor in eq. (7.4) at the first level crossing 



I 


L [fm] 


Q 


qd(f)/dq 


kdSl/dk 





5.34 


0.89 


4.70 


1.12 


2 


6.09 


1.02 


6.93 


-0.09 



physical amplitudes are given by 

Note that Aq and A 2 are real and only their relative sign is observable. Up to this 
sign, the complete information can thus be retrieved from the matrix elements and 
the energy spectrum in finite volume. 

For illustration, let us suppose that the scattering phases Sq are accurately de- 
scribed by the one-loop formulae of chiral perturbation theory [20-22]. The two-pion 
energy spectrum in the subspaces with isospin I and the box sizes L, where the next- 
to-lowest levels in these sectors (the ones with level number n = 1) coincide with the 
kaon mass, can then be calculated. After that the proportionality factor in eq. (7.4) 
is easily evaluated (table 1) and one ends up with 



\A \ = 44.9 x |M |, 
\A 2 \ =48.7x |M 2 |, 
\A /A 2 \ =0.92 x \M /M 2 \. 



(7.5) 
(7.6) 
(7.7) 



As can be seen from these figures, the large difference between the scattering phases 
in the two isospin channels (about 45° at k = k n ) does not lead to a big variation in 
the proportionality factors. In fact, if we set the scattering phases to zero altogether, 
eqs. (4.6)-(4.8) give \Aj\ = 47.7 x \Mj\ for n = 1, which is not far from the results 
quoted above. 

The proportionality factor in eq. (7.4) thus appears to be only weakly dependent 
on the final-state interactions. In particular, if the theory is to reproduce the AI = 
1/2 enhancement, the large factor has to come from the ratio of the finite- volume 
matrix elements Mi. 
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8. Concluding remarks 



Finite-volume techniques have been used in lattice field theory for many years and 
have long proved to be a most effective tool. It may well be that weak transition ma- 
trix elements are also best approached in this way. For two-body decays a concrete 
proposition along this line has been made here, which is conceptually satisfactory 
and which we believe has a fair chance to work out in practice. 

In the case of the physical kaon decays, the proportionality factor relating the 
transition matrix elements in finite and infinite volume turned out to be nearly the 
same in the two isospin channels. This may be surprising at first sight, since the 
interactions of the pions in the isospin state are much stronger than in the isospin 2 
state. One should, however, take into account the fact that the comparison is made 
at box sizes L greater than 5 fm. It is hence quite plausible that the finite-volume 
matrix elements already include most of the final-state interaction effects (such as 
the ones recently discussed in refs. [23-25]). Apart from a purely kinematical factor, 
an only small correction is then required to pass to the matrix elements in infinite 
volume. 

Since the unitarity of the underlying field theory has been essential for our argu- 
mentation, it is not obvious that eq. (4.5) holds in quenched QCD. As usual, however, 
one expects to be safe from the deficits of the quenched approximation when the 
quark masses are not too small and our results should then be applicable. An investi- 
gation of the problem in quenched chiral perturbation theory, following refs. [26,27], 
may be worth while at this point to find out where precisely the unphysical effects 
set in. 

As a final comment we note that the ideas developed in this paper may also be 
applied to baryon decays, such as A — > iV"7r, £ — ► Ntt and S — ► Air, as well as 
to any other decay where the particles in the final state scatter only elastically. 
Depending on the kinematical details, the relation between the finite and infinite 
volume transition matrix elements may, however, assume a slightly different form. 
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Appendix A 



The components of four-vectors in real and euclidean space are labelled by an index 
running from to 3. Bold-face types denote the spatial parts of the corresponding 
four-vectors and scalar products are always taken with euclidean metric, except for 
Lorentz vectors in real space where xy = xo2/o — xv - 

States \p) in infinite volume describing a spinless particle, with mass m and four- 
momentum 



P=(Po,p), Po = V™ 2 + P 2 > °) (A- 1 ) 

are normalized in such a way that 

(p\p') = 2p (2nf5(p-p>). (A.2) 

Particle states in finite volume are always normalized to unity. 

In the centre-of-mass frame, the elastic scattering amplitude of two spinless par- 
ticles of mass m may be expanded in partial waves according to 

oo 

T= WTrWj^i 21 + l)Pi(cos6)ti{k), W = 2y / m 2 + k 2 , (A.3) 

1=0 

where W denotes the total energy of the particles, 6 the scattering angle and Pi{z) 
the Legendre polynomials [28]. Below the inelastic threshold, unitarity implies 

*' = 2Tk ~ ^ ' (A - 4) 

with 5i the (real) scattering phase for angular momentum I. 
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Appendix B 



For all q > the angle 4>{q) is determined through 

3 /2 

tan <f>(q) = - J Q 0(0) = 0, (B.l) 

^ooli;<7 ) 

and the requirement that it depends continuously on q. The zeta function in this 
equation is defined by 

Z 00 (s;q 2 ) = ^= J2(n 2 -q 2 )- s (B.2) 

if Re s > | and elsewhere through analytic continuation. 

Numerical methods to compute the zeta function are described in ref. [17] and a 
table of values of <p(q) is included in ref. [18]. The source code of a set of ANSI C 
programs for these functions can be obtained from the authors. 
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